Theory of De-Pinning of Monolayer Films Adsorbed on a Quartz Crystal 

Microbalance 
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In quartz crystal microbalance (QCM) studies of the friction between an adsorbed monolayer 
film and a metallic substrate, the films are observed to slide relative to the substrate under inertial 
forces of order lQ~^'^dyn per film atom, a force much smaller than all theoretical estimates of the 
force that surface defects are capable of exerting. In this letter we propose, in order to resolve this 
issue, that if the defect potentials have a range of greater than an atomic spacing, the net force on 
a relatively stiff film due to the defects is likely to be extremely small. Line defects (e.g., step and 
facet edges and grain boundaries) as well as more localized defects (e.g., vacancies) are considered. 

PACS numbers: 68.35.Af,62.20.-x 



Quartz crystal microbalance (QCM) studies of mono- 
layers of molecules on metallic substrates [1] provide de- 
tailed information about friction at the atomic level. 
The QCM consists of a quartz oscillator of frequency 
~ IQ^Hz. Monolayer films adsorbed on the metallic sur- 
face of the quartz crystal oscillator are found to slip with 
respect to the surface during a period of the oscillator. 
Small changes observed in the frequency and in the as- 
sociated Q factor allow us to gauge the friction between 
mono-layer and substrate. The amount of dissipation 
generated under most conditions implies that sliding mo- 
tion of more than a lattice constant occurs during a pe- 
riod of the QCM. 

For film atom mass m « lO^^^g, appropriate for xenon 
atoms, QCM frequency uj k, lO^rad/s and amplitude 
A w 100A°, which are appropriate parameters for the 
QCM experiment, the inertial force mw^A is only about 
10~^^ dyn per film atom. Calculations of pinning of a 
monolayer film by defects, based on both perturbation 
theory [2] and on molecular dynamic simulations [3] imply 
that this inertial force per atom due to oscillations of the 
substrate should not be sufficiently strong to overcome 
the pinning resulting from the defects. This is in fact 
what was found in experimental studies by Taborek[4] 
to attempt to reproduce the striking temperature depen- 
dence of the sliding friction found by Dayo and Krim for 
a nitrogen film sliding on a lead substrate, as the tem- 
perature dropped below the superconducting transition 
temperature of the lead [5]. It has been established [6], 
however, that the substrates used in this experiment were 
much more contaminated than those used in Ref. 5. In 
the usual treatment of an elastic medium interacting with 
a disordered potential , the range of the defect poten- 
tial is assumed to be sufficiently small compared to a film 
lattice spacing, so that it can only interact with one film 
atom at a time. In contrast, pinning defects on the sub- 
strate typically extend over one or more atomic lengths. 
For example, the potential produced by line defects (e.g., 
step and facet edges and grain boundaries) clearly extend 
a number of lattice spacings along the length of the de- 



fect, and even a more localized defect such as a vacancy 
extends at least out past its nearest neighbor atoms, be- 
cause the neighboring atoms will displace towards the 
missing atom, resulting in a potential which extends out 
at least that far. It will be argued here that for poten- 
tials, typical of surface defects, that have a long enough 
range to interact with two or more atoms along the slid- 
ing direction, the forces on these atoms tend to cancel 
each other, reducing the force on the film due to a single 
defect to a value much smaller than the maximum force 
that a defect can exert on a single film atom. 

It has been established using simulations |3| that the 
Larkin length 7] , which is a measure of the stiffness of the 
film, is many lattice spacings long implying that the film 
distorts over distances long compared to the mean defect 
potential spacing, and hence certainly long compared to 
the defect potential width. It will now be demonstrated 
that because of the film stiffness, even for defect poten- 
tials with a range of a little more than a film lattice sp 
acing, the forces exerted on the film atoms tend to cancel 
to a great degree. (Incidentally, the simulations reported 
in Ref. 3 used a defect potential of much shorter range.) 
Consider the potential energy of a rigid film of atoms 
interacting with a single defect potential 
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where w(R -I- Ar) is the potential energy of an atom lo- 
cated at the point R -I- Ar in the film due to the defect, 
Ar represents a displacement of the film relative to the 
potential and R is an atomic position in the periodic lat- 
tice of the film. Since this quantity is a periodic function 
of Ar with the periodicity of the film lattice, and hence, it 
can be expressed as a Fourier series 8J in the coordinates 
in the plane parallel to the interface 

^«(R-|-Ar)=^^;(G,z)e^«-^"- (2) 

R G 

where the Fourier coefficient at a distance z from the sur- 
face, v{G, z) = J d^re~**^ ''f (r), where 17 is the unit 
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cell area and G is a reciprocal lattice vector of the film. 
The integral is taken over the whole film. For line defects, 
such as step and facet edges and grain boundaries v(r) 
will vary quite slowly as a function of r along the length of 
the line defect (defined as the distance over which the de- 
fect is relatively straight). As a consequence, we shall see 
that the Fourier coefficient v{G, z) are likely to be quite 
small for general directions of G, such that G has a com- 
ponent along the length of the defect sufficiently large 
compared with 27r divided by the length of the defect, a 
condition which is easy to satisfy for any line defect long 
compared to a film lattice constant. The Fourier coeffi- 
cients will also be very small for defects such as vacancies 
which are localized around a point on the substrate. 

Since the exact form of the defect potential is not 
known, let us illustrate this effect by studying a cou- 
ple of simple smooth potentials which drop off rea- 
sonably rapidly at large distances, in order to develop 
a picture of what one expects for defect potentials 
with a range greater than a lattice spacing. Here we 
will consider primarily a two dimensional Gaussian and 
Lorentzian potential, which represent two extremes, as 
the Gaussian falls to zero very rapidly at large dis- 
tances while the Lorentzian falls off relatively slowly 
(as is evidenced by the fact that its second moment 
is infinite). For an anisotropic Gaussian defect poten- 
tial, v{t) = -yo(^)e-[(^/''i)'+(^/''=)'l, where 61 and 62 
are the range parameters, we obtain, for a triangular 
lattice v{G) = -2yo7r(6i62/(3)i/2a2)e-(G^&?/4+fi?,6^/4) 
For the circularly symmetric two dimensional Loren- 
tian potential v{y) = -Vo{z)[l + {x/hf + (y/&2)^]"\ 
v{G,z) on the right hand side of Eq. (2) is equal 
to -[47r6i62/a'(3)i/2]yoKo(Q), where Q = {Glbj + 
G'^b^Y^'^ and Ko{Q) is the spherical Hankel function 
of the first kind with an imaginary argument [9], 
whose large argument assymptotic form is Kq{Q) k, 
(2/7r(3)^/2e^'3. These forms can be used to model a line 
defect which runs along the x-axis if we take 61 >> 62, 
with hi representing the distance along the defect length 
that one must travel in order to reach a point at which 
the relatively straight section of the defect ends, and 62 
represents the width of the defect potential perpendicu- 
lar to the direction of the line defect. We will assume the 
film's crystallographic axes and sliding direction to be at 
an arbitrary angle with respect to the x and y axes. If we 
wish to model localized defects, such as vacancies or in- 
terstitials, we set 61 &2- The magnitude of the Fourier 
coefficient is a very sensitive function of the range of the 
potential. By taking the negative of the gradient of Eq. 
(2) with respect to Ar, we obtain the force on the film, 
whose maximum value is approximately Gv(G^ z), where 
G is one of the smallest reciprocal lattice vectors. The 
maximum force on the film for the Gaussian defect po- 
tential for &i = 62 = 6 = 1.23a is 0.979 x lQ''^Vfi/b and 
for 61 = 62 = 6 = 0.5a it is 0.0248^0/^- The correspond- 
ing value for the maximum force for the two dimensional 
Lorentzian is 0.00349Vo/6 for 61 = 62 = ^ = 1.23a and 
Qm2Wo/b for bi = b2 = b = 0.5a. These quantities 



should be compared to the maximum possible force that 
these model potentials can exert on a single film atom, 
which is 0.429Vb/^ for the Gaussian and 0.5Vb/& for the 
Lorentzian. We see that for both of these values of b 
the maximum force on the film is much smaller than the 
maximum force that it can exert on a single atom for 
both model potentials and it is a very sensitive function 
of b/a. 

For line defects, the effects can be even more dramatic. 
For example, for a line defect along the x-axis, it is rea- 
sonable to assume that 61 is a number of lattice spacings 
long, and hence, considerably longer than a film lattice 
spacing. Hence, for virtually all orientations of the film 
axes, IGaj^il >> 1, and hence, we can see for the two di- 
mensional Gaussian and Lorentzian potentials discussed 
in the last paragraph, their Fourier coefficients will be 
extremely shall, and hence we see from Eq. (2), the de- 
pendence of the interaction of the potential with the film 
with the defect on Ar will be extremely small, implying 
extremely small forces on the film. 

Other model defect potentials give similar results. For 
example, the potentials —Vo/[cosh{x/bi)cosh{y/b2)] and 
— Vo/[(l + (•^;/fei)^"(l + (y/^2)^"'], where n is any integer, 
also resiilt in Fourier coefficients which fall off as expo- 
nential functions for reasonably large values of |Ga; |6i and 
\Gy\b2- The latter potential with a sufficiently large value 
of n is a good model for a line defect whose potential is 
nearly constant along most of its length and only has sig- 
nificant variation near its ends. The only requirement in 
order to get a Fourier coefficients which fall off rapidly at 
large values of its argument (i.e., exponentially or better) 
is that the potential due to the defect be smooth. 

At any temperature, there will be lattice vibrations. At 
a given instant of time, a lattice vibration of wavevector 
q will add a term Acos(q- (r-l-R) —ojt) to R in Eqs. (1), 
where A is the amplitude (assumed to be much smaller 
than a lattice constant a) and w is the frequency of the 
vibrational mode and t is the time. Then, expanding v 
in a Fourier series, Eq. (1) becomes 

J2 J2 z)e*-(^+^'-)e*-^'=°^('»-('-+^)-'^*), (3) 

r k 

where N is the number of atoms in the film. Expand- 
ing the second exponential in Bessel functions Eq. (3) 
becomes 

00 

iV-i Y. ^""M^ ■ A)e*(''-"i)-(^+^"')e*"'^*, 

R k n=a 

(4) 

which when summed over R becomes 
^ ^ z" J„((G - nq) • A)v{G - nq, 2)e'(G-Ar-na;t) _ 

G n 

For the n=l term, which should be a good approximation 
for the small values of | (G — nq) • A| characteristic of lat- 
tice vibrations, we find that v{G — q, z) for the Gaussian 
potential now contains a term e~l^~''l ^ Z^, which can be 
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considerably larger than e~'^ ^ which occurs for n=0, 
by a sufficient amount to make the former term dominate 
over the latter for values of q which are comparable in 
magnitude to the smallest values of G, even though it is 
multiplied by J„((G — q) • A), which is much less than 1. 
As the former term is proportional to A which oscillates 
in time for a lattice vibration, this term produces an os- 
cillating term in the force on the film, which can be larger 
in magnitude than the force on the film that would occur 
if there were no vibrations [i.e., the n=0 term in Eq. (5)]. 
Thus at nonzero temperatures, the lattice vibrations can 
produce oscillating forces acting on the part of the film in 
the range of the potential, which dominate over the static 
pinning force on the film due to the potential under con- 
sideration, implying that at nonzero temperatures, the 
film will not be pinned by the defects because the film is 
pushed out of its total potential minimum by this force. 

In the previous paragraphs, the film was taken to be 
completely rigid, and the distortion of the film result- 
ing from its interaction with the defect potential was ne- 
glected. In order to examine the validity of this approx- 
imation, let us now consider the small distortion of the 
film resulting from the potential, in order to determine 
whether it will destroy the reduction of the potential's 
interaction with the film discussed in the previous two 
paragraphs. The displacement of a film atom located at 
point R is given in lowest order perturbation theory by 

u(R) « -^G(R-R') ■ V'u(R',^), (6) 

R' 

where G(R — R') is the elasticity Green's function (8(|. It 
is shown in Ref. @ that 

£ - „jk-(R-R') 

G(R - R') = N-^ Z^iv^ ' 

where Wq, (k) is the frequency of the a*'' vibrational mode 
of wavevector k and ek.a is the unit vector which gives its 
polarization. Substituting Eq. (6) in Eq. (6), we obtain 

u(R) « -N^\ ''"''vtr'' ■ 

If we add u to R in Eq. (1), linearize in u and substitute 
for u from Eq. (9), we find that the first order term in 
an expansion in powers of u for the force on the film is 
given by 

y [(G-k).g, ][k..V„] _ i^^-(^)^.k.Ar^ 

(9) 

We see from Eq. (9) that since we found earlier that 
v{\i, z) becomes negligibly small rapidly once k becomes 
larger than 1/b, only Fourier components of wavevector 
larger than 1/b will make a significant contribution to 
Eq. (9). We shall now illustrate that the modulation of 
the film resulting from the defect potential should have 



0.000040 




-0.000040 I ' ' ' ' ' ' ' 1 

-4 -2 0^2 4 

s(a) 

FIG. 1: The total force due to a Gaussian potential acting on 
a rigid lattice of atoms of spacing a is plotted as a function 
of the displacement of the film s in the lower amplitude curve 
for b=l. 236067978a. The total force acting on a modulated 
lattice is shown in the higher amplitude curve. 



little effect on our conclusion that the defect interacts ex- 
tremely weakly with the adsorbed film with a simple cal- 
culation. Consider a stiff two dimensional square lattice 
of atoms of lattice constant a interacting with a Gaus- 
sian potential —Voe~^ , were r is the distance from the 
center of the potential. We choose the range parameter 
of the potential, b=1.23a. For this choice of parameters 
the potential range is about one and a half lattice spac- 
ings. The lattice is slid a distance s and the total force 
due to this potential acting on it is plotted as a func- 
tion of s (the lower amplitude curve in Fig. 1). The 
position of an atom in the lattice is then displaced by 
an amount 0.01acos((27r/A)a;„) in the x-direction, where 
the position of a film atom is {xn,ym) = {na,ma) where 
n and m are integers, to simulate the effect of this static 
modulation of the film. There is no displacement Aa; in- 
cluded because a modulation produced by the potential 
does not slide relative to the potential. Here, A is cho- 
sen to be 4b (in the range of the shortest A value which 
makes a significant contribution in Eq. (10)), so that half 
of a wavelength extends across the width of the potential 
well. The results are given in Fig. 1. As can be seen the 
force on the film when such a modulation is present is 
comparable to that in the absense of a modulation, and 
is much less than the maximum force that a defect can 
exert on a film atom. 

We have considered the interaction of a single defect 
with the film and have demonstrated that this can easily 
be much weaker than the maximum possible interaction 
of that defect with a single film atom. When we consider 
the interaction of the film with the random distribution 
of defects over the substrate that is typically found, we 
must consider the fact that in the weak pinning regime 
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(in which the film is not able to distort enough to mini- 
mize its interaction with all of the randomly distributed 
defects) , the forces due to the defects do not act in phase 
0. As a consequence, the net force per atom on the film 
gets reduced in addition by a factor of the square root of 
the number of defects within a Larkin domain, which is 
proportional to the ratio of the Larkin length and a film 
lattice constant. Since for a two dimensional solid the 
Larkin length is equal to the ratio of the elastic constant 
and the interaction with a sing le defect [HI, 

the over- all 

interaction of the defects with the film is proportional to 
the square of the defect interaction, rather than being a 
linear function of it. Furthermore, our single defect in- 
teraction estimates are in fact an upper bound on this 
quantity. 

In conclusion, we have proposed a mechanism which 
may explain why local defects, which must certainly be 
present on even the smoothest surfaces, might not pre- 
vent a stiff monolayer film from sliding under the ex- 



tremely weak inertial forces that occur in a quartz mi- 
crobalance experiment, if the defect potential has a range 
greater than about 1.5 film lattice spacings. Our pro- 
posed mechanism is supported by a simple model calcu- 
lation, based on the fact that each defect can interact 
with several film atoms at a time. As a consequence, 
the resultant forces on the film are much less than the 
maximum force provided by the defect potential. Lattice 
vibrations are able to provide an oscillating force which 
is larger than the net static force due to the defect po- 
tential acting on the film. This can lead to a high degree 
of thermal activation of the sliding of the film. 
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